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ABSTRACT 

Recent deep Hubble Space Telescope WFC3 imaging suggests that a majority of compact 
quiescent massive galaxies at z ~ 2 may contain disks. To investigate this claim, we have 
compared the ellipticity distribution of 31 carefully selected high-redshift massive quiescent 
compact galaxies to a set of mass-selected ellipticity and Sersic index distributions obtained 
from 2D structural fits to ~40, 000 nearby galaxies from the Sloan Digital Sky Survey. A 
Kolmogorov-Smirnov test shows that the distribution of ellipticities for the high-redshift galaxies 
is consistent with the ellipticity distribution of a similarly chosen sample of massive early-type 
galaxies. However the distribution of Sersic indices for the high-redshift sample is inconsistent 
with that of local early-type galaxies, and instead resembles that of local disk-dominated 
populations. The mismatch between the properties of high-redshift compact galaxies and 
those of both local early-type and disk-dominated systems leads us to conclude that the basic 
structures of high-redshift compact galaxies probably do not closely resemble those of any single 
local galaxy population. Any galaxy population analog to the high-redshift compact galaxies 
that exists at the current epoch is either a mix of different types of galaxies, or possibly a unique 
class of objects on their own. 



Subject headings: galaxies: evolution — galaxies: formation — galaxies: fundamental parameters 
galaxies: statistics — galaxies: structure 



1. Introduction 

With the discovery of very compact mas- 
sive quiescent galaxies at redshift s z ~ 2 (e.g 
Daddi et all 120051: iToft et alJliool iTruiillo et al 
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2008 
2009; 



van Dokkum et all 12008 
20101 : 



Damianov et al 



Saracco et all 2011; 



, Williams ct al 
ISzomoru et all 2011 . and references therein), 
galaxy formation and evolution models have 
been challenged to explain the increase in size 
of these systems with decreasing redshift and 
the fact that similarly compact massive galax- 
ies are almost completely absent from the lo- 
cal Universe (ITruiillo et all 120091 : ITavlor et all 
2010l : IValentinuzzi et al.l l2010ll Compared to 



present-day galaxies of similar (stellar) mass 
(M* ~ 10 11 Af Q ), these high-z compact galaxies 
or 'red nuggets' are a factor of ~ 5 smaller (e.g., 
iDamjanov et al.ll2011 1 with half-light or effective 
radii R e < 1 kpc, while their stellar mass densities 
within the effective radius are on average an order 
of magnitude higher (although this difference is 
smaller for the stellar mas s densities with i n the 
central kiloparsec, see e.g. IBezanson et al. 2009; 
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Saracco et al.l 120121) . Based on the low fraction 
of close pairs among quiescent galaxies observed 
at < z < 2 (Bell et al.ll2006t iBundv et al.ll20Qg : 



De Propris et all l2010i: IWilliams fc Quadril [2011 



Man et al. 20121) and the small number of near 



equal mass mergers produced in N-body simu- 



lations (e.g., Shankar et al. 2010l ). it seems that 



major mergers can only be partly responsible for 
the observed size growth. Additional s ecular pro- 
cesse s , such as adiabatic e xpansion (IFan et al 



20ld : Damianov et al. 20091) and / or a series of 
mino r me rgers (jNaab et all 120091 : iHopk ins et al.l 
20101 : lOser et a.1.1 l2012h . might be needed to ex- 
pand these compact systems. 

Recent deep, high-resolution images taken with 
the Wide Field Camera (WFC3) and NIC2 cam- 
era on board the Hubble Space Telescope (HST) 
have shown that based on observed elliptiti- 
cies and Sersic profile fits, there may be indi- 
cations that many of th ese compact high-redshift 
galaxies contain disks ( van Dokkum et all 2008: 



van der Wei et al .2011). Indeed, Ivan der Wei et al 



(|201ll ) claim based on their sample of 14 quiescent 
massive compact galaxies, that the majority of 
these systems are disk dominated. Based on the 
structur al parameters de r ived f rom ground-based 
imaging, IWhitaker et al.l (|201 2| ) have recently re- 
ported a small decrease in the average axial ratio 
with redshift (i.e., mildly increased prominence of 
the disk component) for a sample of post-starburst 
galaxies found at < z < 2. These findings would 
have interesting consequences for the possible for- 
matio n and evolution scen arios of these galaxies 
Weinzirl et al.l 120111) . I n the present pa- 



(e.g., 

per w e investigate the claim in Ivan der Wei et al 



( 201ll ) by comparing the ellipticity and Sersic dis- 
tributions of these objects to the corresponding 
distributions for nearby massive disk-dominated 
galaxies, as well as massive bulge-dominated early- 
type galaxies. To defi ne our high-redshif t sam - 
ple we complem ent the |van der Wei et all ( 201lh 
sample with the iDamianov et al] (|201ll ) synthesis 
of published structural data for high-z compact 
galaxies with secure spectroscopic redshifts. At 
low redshift, we characterize the ellipticity and 
Sersic index distributions of nearby galaxies using 
the local morphology sample recently published 
bv lSimard etaD (1201 ll) . 



2. Data 

Our sample of high-redshift massive compact 
galaxies consists of 14 galaxies with photometric 



redshi fts z > 1.5 taken from Ivan der Wei et al 
([201 ll ). augmented by object s selected from 



the c ompendium published by IDamianov et al 
(|201ll ). All 69 galaxies in the Damjanov et al 



sample have spectroscopic redshifts, HST-based 
images and published ellipticities and Sersic in- 
dices. Around 1/4 of these galaxies are at z > 1.5, 
have stellar masses M* > 10 la8 M Q , and are com- 
pact (R e < 2 kpc) and quiescent (not actively 
forming stars). These 17 objects are added to our 
high-redshift sample, resulting in a total sample 
size of 31 galaxies. Galaxies in the Damjanov 
et al. sample span a mass range from 10 9 ' 8 Mq to 
lO^^Mcn, with a median of 10 u l M ra . Galaxies in 



the Ivan der Wei et al.1 (|201lt ) sample all have stel- 
lar masses larger than 10 11 ' 1 M Q , with a median of 
10 113 M Q when scaled to a Salpeter (1955) initial 
mass function to be consistent with the Damjanov 
et al. sample. The two samples therefore include 
objects with similar stellar masses and in the same 
redshift range, and all have R e < 2 kpc. 

Our local galaxy sample is selected from 
the morphological c atalog recently published by 



Simard et al 



( 201 lh . This catalog is based on 
two-dimensional structural fits (Sersic bulge+ ex- 
ponential disk decompositions, as well as single 
Sersic fits) in the g and r bandpasses for 1.12 mil- 
lion local galaxies taken from the Sloan Digital 
Sky Survey. We restrict our comparison to galax- 
ies with redshifts z < 0.1 in order to minimize 
the impact of seeing on the obs erved shapes. As 
described in Simard et al. (2009), a robust sample 
of early-type galaxies can be selected from this 
parent catalog by requiring that galaxies have a 
bulge fraction B/T > 0.5 and an image smooth- 
ness 52 < 0.07£1 Imposing these cuts results 
in a sample of ~ 127,000 early-type galaxies, 
which defines our local early-type galaxy parent 
sample. Our local disk-dominated sample was 
chosen solely based on bulge-to-total ratio, with 
B/T < 0.5. Stellar masses, M„, were estimated 



using the approximation presented by iBell et al 



1 As described below, we experimented with changing these 
parameters to select even cleaner subsets of early-type 
galaxies. All of the results in the present paper are robust 
to modest changes in these parameters. 
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(|2003h : log(M*/M ) = -0.499+1. 519(M S -M r 



0.4 x (M g - Af 90 ) - log(0.7) where M g and M r 
are rest-frame absolute magnitudes in g-band and 
r-band, respectively, and M 90 — 5.12 is the g- 
band absolute magnitude of the sun. A ISalpeter 
(Il955l) initial mass function (IMF) was assumed in 
deriving this relation. The relative masses deter- 
mined using this formula are likely to be accurate 
to 0. 1 dex over the re dshift range of our local sam- 
ple ( Bell et al.l 120031 ). Fair comparisons between 
the high-redshift and low-redshift samples require 
careful matching of the mass ranges of the sam- 
ples, and most of the comparisons in this paper 
will be based on mass-selected subsets of our local 
galaxy sample. Note that our selection of our lo- 
cal samples is based on morphology only, and does 
not put any restrictions on star formation activity 
or quiescence. 



3. Analysis 

Figure [1] presents the apparent ellipticity dis- 
tributions for our high-redshift massive compact 
galaxy sample. The ellipticity is as usual defined 
as e = 1 — b/a, with b/a the observed axis ra- 
tio of the images, such that round galaxies have 
e = 0. These ellipticities (and also the Sersic in- 
dices that we introduce below) have been mea- 
sured on the total galaxy image using the two- 
dimensional fitting routin e GIM2D, so prio r to 
bulge-disk decomposition (jSimard et al.ll201ll ). In 
the left-hand panel, ellipticity distributions for 
our local early-type galaxy samples have been su- 
perposed on the high-redshift distribution. The 
local early-type galaxies samples consist of ob- 
jects with stellar masses > 1O 1O8 M (solid red 
line) and early-type galaxies with stellar masses 
> 10 n M Q (dotted purple line). The > 10 11 M 
curve provides the most meaningful comparison 
to the compact massive galaxies sample, but both 
local galaxy distributions are shown in order to il- 
lustrate the robustness of our conclusions to small 
(0.2 dex) variations in the mass threshold. A 
Kolmogorov-Smirnov (K-S) test was used to deter- 
mine whether the low-redshift distributions differ 
significantly from the high-redshift distribution, 
and the corresponding p-values from the test are 
shown as labels in the figure. The high-redshift 
and local early-type galaxy ellipticity distributions 
differ at a less than la level and are therefore 
statistically indistinguishable. Our conclusion re- 



mains the same if we choose a more conserva- 
tive cut (B/T > 0.7) in defining our local early- 
type galaxy sample (p = 77% and p = 67% for 
M* > 1O 1O - 8 M and M* > 10 U M Q respectively). 

Although the high and low-redshift ellipticity 
distributions presented in the left panel of Figure[T] 
appear consistent, the small number of galaxies in- 
cluded in the high-z sample raises the question 
whether there are any distributions at all that 
could be ruled out using such a small sample. It is 
thus rather interesting to show that one distribu- 
tion that can clearly be eliminated is that of local 
disk-dominated galaxies. The right-hand plot in 
Figure [T] shows the results of K-S tests compar- 
ing our sample of 31 high-redshift massive com- 
pact galaxies to the ellipticity distributions of sim- 
ilarly massive disk-dominated (B/T < 0.5) galax- 
ies from our local sample. The ellipticity dis- 
tribution for this local population decreases for 
e < 0.2, indicative of earlier findings that lo- 



cal d isks are not intrinsically circular (e.g. iRvden 
20041 ). There is only a 5% probability that the 



high and low-redshift ellipticity distributions are 
drawn from the same population. For a more con- 
servative cut (B/T < 0.35), this probability drops 
to 2%. Therefore, on the basis of ellipticity dis- 
tributions alone, one is driven to conclude that 
the population of compact objects is dominated 
by early-type galaxies. Our conclusion remains 
unchanged if we select from our local samples only 
those galaxies that have effective widths (EW) of 
the [Oil] line less than 5A, indicating a quies- 
cent (non-star forming) nature. For local mas- 
sive (M* > 10 11 Mq) quiescent early-type galax- 
ies, p = 82%, while for the local massive quiescent 
disk galaxies, the probability that they come from 
the same distribution as the high-redshift com- 
pacts stays low at p = 11%. However, is the same 
conclusion reached when other structural data are 
also considered? 

A simple structural parameter that nicely com- 
plements ellipticity is the Sersic index, n, which 
is essentially a measure of galaxy compactness. 
This parameter is strongly correlated with Hubble 
type, and nearby systems with n < 2 are generally 
disks, although low-mass early-type galaxies and 
dwarf sph eroidals also mostly have l ow Sersic in- 
dices fe.g. iGraham Sz Guzman! 12003 ) . Repeating 
the K-S tests already described, but now using his- 
tograms of Sersic index instead of ellipticity, pro- 
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Fig. 1. — (Left:) The ellipticity distribution for 14 compact massive galaxies at z ~ 2 from van der Wei 
et al. (2011), complemented by 17 high-redshift massive compact galaxies from Damjanov et al. (2011) 
(broad grey bins). These data are compared to the distributions of local massive early- type galaxies at 
z < 0.1 from Simard et al. (2011). Solid red and dotted black lines show the ellipticity distributions of local 
early-type galaxies with > 10 11 Mg and M* > 10 10 8 M Q , respectivly. The labeled percentages are the 
corresponding p-values from a Kolmogorov-Smirnov test which gives the probability that the high-redshift 
data and low-redshift data originate from the same distribution. (Right:) As for the panel at left, but now 
comparing the ellipticitity distribution of the high-redshift massive compact galaxies from van der Wei et al. 
(2011) and Damjanov et al. (2011) to the distribution of massive nearby disk-dominated galaxies at z < 0.1. 



vides us with a different picture of the structural 
composition of the high-redshift compact popula- 
tion. As shown in Figure 2, the distribution of 
Sersic indices of local early-type galaxies is clearly 
inconsistent with that of high-redshift compact ob- 
jects (see left panel Figure [5]). The p- value re- 
turned by the K-S test is ~ 10~ 15 , reinforcing 
the visual impression that the probability of these 
samples being drawn from the same underlying 
distribution is negligible. Interestingly, the K-S 
test does return a much larger p- value (60%) when 
we compare the Sersic index distribution of the 
local massive disk-dominated galaxy sample with 
the high-redshift one, though only for the most 
massive (M* > 10 Mq) systems (right panel of 
Figure [5]). Selecting only quiescent (EW[OII] < 
5A) local galaxies does again not alter our con- 
clusions: for local massive quiescent early-types 
p ~ 10~ 14 , although for local massive disk galax- 
ies p = 33%. 

Sersic indices are challenging to measure in 
high-redshift galaxies, as possible outer wings of 
surface brightness profiles may be too faint to ob- 



serve, biasing the observed Sersic indices towards 
lower values. As a consequence, the systematic 
errors on these measurements could be significant 
and influence our conclusions above. We therefore 
investigated the effect that measurement errors in 
the Sersic index could have on our results with 
Monte Carlo simulations. The errors in the Dam- 
janov et al. sample were determined with Monte 
Carlo simulations which incorporated all the sys- 
tematic and random er rors that they were ab le to 
identify (for details see lDamianov et"al"1l2009l), ex 



changes in background (Ivan Dokkum et al 



or GALFIT statistical errors (Dad di et al 



2008) 
20051) . 



We do n ot have errors for the Ser sic indices pre- 
sented in Ivan der Wei et al] (|201 lh , but adopt for 
these galaxies the median value of the errors in 
the Damjanov et al. sample, which is 0.5. We 
run the K-S test 10,000 times on our sample, each 
time adding Gaussian noise to our observed Sersic 
measurements based on these errors. 

For the K-S tests comparing the Sersic indices 
of the high-redshift sample with the highest mass 
(M* > 10 11 Af©) local early-type galaxies, we ob- 



4 




Fig. 2. — As Figure HJ but now comparing Sersic index distributions. The excess of objects with n > 7 
consists mostly of galaxies with a nuclear source or a bar, or a combination thereof (Simard et al. 2011). 



tain p ~ 10 -12 . Even if we place all Sersic indices 
at their one-sigma upperlimit, we find p ~ 10 . 
We therefore conclude that the missmatch be- 
tween Sersic indices of our compact high-redshift 
and local early-type sample is not a result of possi- 
ble underestimation of the measured Sersic indices 
at high redshift. 

4. Discussion 

Compact massive quiescent galaxies at z ~ 2 
have an ellipticity distribution consistent with 
that of local massive early-type galaxies, but their 
Sersic index distribution seems to rule out the 
possibility these high-z objects are all early-type 
galaxies. The Sersic index distribution of the high- 
redshift galaxies is fairly similar to that of local 
disk populations, though the ellipticity distribu- 
tions are very different. What are we to make of 
this conundrum? One possibility is that the pop- 
ulation of high-redshift compact massive galaxies 
is a mixture of bulge and disk dominated galax- 
ies. Such a model might be made to work, but in- 
spection of Figures 1 and 2 suggests that it would 
not be very straightforward to mix ellipticity and 
Sersic indices for early-type and disk-galaxies in 
such a way as to have bulge-dominated galax- 
ies dominate one property, while disk-dominated 
galaxies dominate the other. Another possible 
problem with this model is that it would require 
galaxies, that have different structural properties, 



to all change in size with cosmic epoch, which 
seems hard to credit unless the growth mechanism 
is somehow extrinsic to the gal axies themselves 



a hypothesis briefly discussed in lDamianov et al 



20111 ). 



An alternative, and in some ways more natu- 
ral explanation for the discrepancy between com- 
bined structural properties of our samples at z ~ 2 
and z ~ might simply be that the high-redshift 
galaxy population has unique properties, that are 
not present in the local galaxy population because 
local galaxies are poor analogs to the high-z sys- 
tems. In other words, the high-redshift compact 
galaxies may truly be a class on their own. 

In Figure |3] we explore these alternatives a little 
further by plotting ellipticity versus Sersic index 
for our high-redshift sample on top of contours 
showing the distribution of these same parame- 
ters in local galaxies. Different panels in this figure 
correspond to different mass bins. Table [T] summa- 
rizes the results for the various K-S tests that we 
ran on all these samples. None of the local galaxy 
samples provide a good fit to the observed distri- 
butions at high redshift, and the fact that also the 
lower-mass early-type galaxy samples fail to do so 
rules out the possibility that our results are due 
to errors in assigning stellar masses to the galax- 
ies (perhaps as the result of an evolving IMF). 
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Fig. 3. — Comparing elliptiticity and Sersic index distributions from local early-type galaxies (contours) to 
those of the high-redshift massive compact galaxies of van der Wei et al. (20 f f , black filled triangles) and 
Damjanov et al. (2011, black stars). Contours are normalized and smoothed, increasing in number density 
from yellow to red in logarithmic steps. Different panels show different mass ranges for the local early-type 
galaxy sample, denoted in the upper right corner or each panel. The last panel (lower right) shows the 
distribution for local massive disk- dominated galaxies. 



Interestingly, a simple likelihood analysifd shows 
that the disk-dominated sample (most-right bot- 
tom panel) has a higher likelihood than the early- 
type local samples, when taking both the ellip- 
ticity and Sersic distributions into account simul- 
taneously. This likelihood analysis however only 
tells us that one distribution is a better fit than 
the other, but not that this fit is necessarily a good 
one. 

Figure [3] highlights another important point. 
Although no local bivariate ellipticity- Sersic in- 
dex distribution is a good fit to the high-z 
data, there certainly does exist a strong over- 
lap in parameter space between the high-z pop- 
ulation and the local sample. This fact, to- 
gether with inspection of the images presented 



2 We calculate the likelihood L as logL = 

P(ej, nj)dedn, for our sample of N = 31 high- 

z galaxies, with P(e,n) defined by the normalized 
distribution of the local galaxy sample. 



in van der Wei et al. I(l201lh. btt also lDaddi et al 
( 20051) :lLonghetti et all (12001 : iToft et al l (| 2007 h : 
van Dokkum et al.l (120081): ICimatti et al.l (120081 ) 



Damianov et al. ( 2009t h iRvan et al ( 2010h : ISzomoru et al 
(|201ll ). makes it seem fairly clear that some of 
these objects a r e ind eed disks. Recent results of 



Weinzirl et al.l (|2011l ) shown that although the 
majority of all massive galaxies at 2 < z < 3 show 
possible disk-like str uctures with n < 2 (see also 



Buitrago et all 1201 lh , most (~ 72%) of these ob- 
jects are extended {R e > 2 kpc). On the other 
hand, half of the population of compact (pre- 
dominantly quiescent) objects do not show any 
disk-like features, and for the highest mass bin 
(M* > 1O U M ) only - 38% has n < 2. 

Our results lead to a prediction that can be 
used to test the formation mechanism(s) of mas- 
sive galaxies. Major gas-rich mergers are be- 
lieved to produce compact spheroidal remnants 
with light pr ofiles characterized by high Sersic 
indices (e.g., Wuvts et al. 20101 : Bournaud et al 
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20111 ). Additional gradual and slow gas accre- 
tion is required to b u ild up a disk componen t 



(jWeinzirl et al.l l201ll : Ivan der Wei et al.l I2OIII ). 



One intriguing possibility is that compact qui- 
escent galaxies at high-redshift consist of thick 
disks, possibly for med by star-forming turbulent 
disks at z > 2 (e.g. Elmegreen fc Elmegreenl 2006L 



Bournaud et alJl2009l ). The thick disk would then 
be responsible for the mismatch with the ellipticity 
distributions of the local disk-dominated popula- 
tion, with the disk thickness driving the lack of 
objects with extreme ellipticities (e > 0.7) com- 
pared to modern spirals, whose light is dominated 
by thin disks. 

It is worth emphasizing that if quiescent pop- 
ulations at z = and z ~ 2 have very different 
galaxy light profiles, which is one possible interpre- 
tation of our data, then the subsequent evolution 
of compact quiescent galaxies bearing disks would 
need to involve a drastic change in their Sersic in- 
dices. In a scenario in which their size evolution 
is driven mainly by a series of minor mergers, the 
Sersic indices of compact disk galaxies can change 
by adding a small amount of l ow density stel- 



lar mass to their outskirts (e.g., iBournaud et al 
2007). At this point it is not clear whether that 
change would be dramatic enough to turn com- 
pact disks int o systems resemblin g the spheroids 
found locally (jWeinzirl et al.ll2011l ). 



5. Conclusion 

We have compared the observed ellipticity dis- 
tribution of massive, compact, quiescent galaxies 
at high redshift (1.5 < z < 2.5) to those of lo- 
cal [z < 0.1) early- type galaxies, and conclude 
that the two distributions are statistically indis- 
tinguishable. In addition, we show that the ellip- 
ticity distribution of our high-z sample is inconsis- 
tent with that of a local massive disk-dominated 
sample. On the other hand, the Sersic index distri- 
butions of compact high-z galaxies and local early- 
type galaxies are not consistent, and the Sersic in- 
dex distribution of local disk galaxy samples pro- 
vides a better match to the high-z data. We con- 
clude that either (a) high-z galaxies are a compos- 
ite population of disks and bulges (which presents 
the troublesome possibility that galaxy sizes are 
growing for all structural types); or else, (b) that 
the shapes and light profiles of high-redshift mas- 



sive compact galaxies are unlike those of any local 
galaxy sample, and their structure changes over 
time. In the latter case, these objects constitute a 
new class of galaxies on their own. 

We thank the referee for his/her constructive 
comments, which helped to improve the presenta- 
tion of this paper. This work makes use of the 
Sloan Digital Sky Survey (www.sdss.org). 
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Type 


Mass Range 


N 


ellipticity 


Sersic index 




(log Mq) 




p- value 


p- value 
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11.0 - 
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0.82 


0.00 




10.8 - 
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0.00 
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10.2 - 10.4 
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0.64 


0.00 


disk 


11.0 - 


30,098 


0.05 


0.60 




10.8 - 


56,763 


0.03 


0.07 



Table 1: Sample sizes and Kolmogorov-Smirnov test p- values (probability that two samples were drawn 
from the same underlying distribution) for comparisons in ellipticity and Sersic index between the local and 
high-redshift sample, for several mass bins. The local early- type sample is chosen such that B/T > 0.5 and 
smoothness parameter S2 < 0.075; the local disk sample is selected such that B/T < 0.5, with no additional 
constraints on smoothness. See text (§2) for more details. 
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